Parameter Estimation with Mixed-State Quantum Computation 
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We present a quantum algorithm to estimate parameters at the quantum metrology limit using deterministic 
quantum computation with one bit. When the interactions occurring in a quantum system are described by a 
Hamiltonian H = 6 Ho, we estimate 9 by zooming in on previous estimations and by implementing an adaptive 
Bayesian procedure. The final result of the algorithm is an updated estimation of 8 whose variance has been 
decreased in proportion to the time of evolution under H. For the problem of estimating several parameters, we 
implement dynamical-decoupling techniques and use the results of single parameter estimation. The cases of 
discrete-time evolution and reference frame alignment are also studied within the adaptive approach. 
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I. INTRODUCTION 

Quantum mechanics provides new resources that allow us 
to determine physical properties at the highest possible accu- 
racy established by generalized uncertainty relations 01 0. HI] . 
Exploiting quantum coherence enables us to estimate parame- 
ters 0, Et] and expectation values of observables J(J with bet- 
ter resource scaling than classically possible. In this paper, we 
are interested in the estimation of interaction parameters (e.g., 
external fields), when the interaction acts independently on n 
quantum subsystems in a probe 0). We quantify our resource 
of interest by N = nT, given by the product of the number 
of subsystems and the interaction time T. The standard quan- 
tum limit (SQL) precision in the estimation of an interaction 
parameter is of order 0(l/y/~N), achievable with O(N) inde- 
pendent measurements at a fixed T. The optimal precision for 
such an estimation, however, is given by the Heisenberg limit 
and is known to be of order 0(1/ N). Achieving it requires 
the preparation of entangled quantum state in the probe ||4[] ■ 

We are interested in estimating parameters at the so called 
quantum metrology limit (QML). This can be obtained by a 
series of estimations performed at different interaction times, 
while keeping the size of the probe fixed JH |^] . If for a total 
interaction time T, the precision of the estimation is of or- 
der 0(1/T), we say that the QML has been achieved. This 
sequential protocol, which is the one exploited in this paper, 
does not require quantum entanglement in the input state, al- 
though the response to uncorrected decoherence, for an un- 
constrained interaction time, is the same as that of the entan- 
gled protocol @]. Any method that allows us to achieve the 
QML provides clearly an improvement over the SQL, since 
for the same amount of resources (i.e., N = nT), the returned 
precision can be highly enhanced. 

Quantum methods (algorithms) designed to beat the SQL 



could have a wide range of applications, from highly sensi- 
tive magnetometry ifioll to atomic clock synchronization Ultl . 
In addition, phase estimation, a problem related to parameter 
estimation, is one of the cornerstones of Quantum Computa- 
tion lfl2i [l3ll . In this paper we show that the QML can be 
achieved in some cases even if the initial state is the com- 
pletely mixed state of all except one of the quantum systems, 
avoiding the complexity associated with initial pure, entan- 
gled, state preparation. Although here we consider multi- 
qubit probes, generalization of our algorithms to higher di- 
mensional systems is straightforward. 

Specifically, we use deterministic quantum computation 
with one bit (DQC1), which was initially described in 
Ref. Hill in the context of high temperature ensemble quan- 
tum computation using liquid-state NMR techniques lfl5ll . Al- 
though less powerful than the standard model of quantum 
computation, DQC1 is believed to outperform the classical 
probabilistic computational model flrjl . In DQC1, the initial 
state po of a set of n + 1 qubits corresponds to having the first 
(ancilla) qubit a in the pure state |0) a , while the state of the 
remaining n qubits (probe) is completely mixed. That is, 
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The state po is then unitarily evolved and DQC1 returns a 
noisy expectation value of a Pauli operator on the ancilla 
qubit. If the evolution is performed by applying a unitary op- 
eration controlled by a (i.e., a controlled-?/ or C U operation), 
DQC1 allows us to estimate the renormalized trace of U at a 
certain, fixed precision (Fig. [T}. 

We assume that a single run of the DQC1 algorithm returns 
an unbiased renormalized trace estimation with normal dis- 
tribution JV(tr[f7]/2™ A 2 ) of standard deviation A fl7il. Of 
course, A can be reduced by a factor of \/~K (i.e., SQL) if the 
algorithm of Fig. [I] is repeated K times. In fact, this is the 
situation in NMR where repetition reduces the signal-to-noise 
ratio (SNR) at the SQL. 

Consider now, for example, the typical case where an un- 
known external magnetic field interacts with the n qubits of 
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•ioS)=tr[Cn/2» 



FIG. 1: DQC1 circuit for estimating the trace of a unitary U. The 
filled circle denotes that U acts on the probe when the state of the 
ancilla (control qubit) is |1) , and H is the Hadamard gate. (cr%) and 
(<jy) are the expectation values of the corresponding Pauli operators 
on the ancilla a. 



the probe (Fig. [T}-. determining an interaction Hamiltonian of 
the form 



(2) 



We seek to estimate 9. When the probe interacts with the field 
for time T, the n-qubit state is evolved by applying the corre- 
sponding (unitary) evolution operator W'(T) 
placing C U by C W'(T) in Fig.Q] the final (n + 
right before the measurement is 



-iff T 



Re- 

1)— qubit state 



Pi = ^pr[l )a( l ® fl « + |!>a a <0| ® W(T) + (3) 
|0) a a (l|®^'t(T) + |l) a a (l|®ll„]. 

Using the trace properties of the Pauli operators, we obtain 
(o%) = tr \pfO%] = (cos(0T)) n . If A x > denotes the stan- 
dard deviation in the estimation of (&%), a first approximation 
error formula determines 



A, 



(4) 



with Ag the uncertainty in the estimation of 9. 

Let |+) denote the single qubit state (|0) + |1))/V2. 
The output signal of the previous algorithm is 
|(+ 1 ...+„ /2 |vr(T)|+ 1 ...+ n/2 )| 2 , which is the 
probability of measuring |+i-..+ n /2) after its evolution 
under W'(T). Thus, the output precision of Eq. (0]i is upper 
bounded by Ag = 0[1 / (y/nT)], and the Heisenberg limit is 
not achieved when scaling n. Nevertheless, for fixed (small) 
n, we obtain Ag = 0{1/T), yielding the QML. 

The previous estimation method has some important dis- 
advantages. The first one concerns the use of the controlled 
C W'(T) operation which, due to technological difficulties or 
to the nature of the problem, may be impossible. In fact, 
this is the case in reference frame alignment, as we discuss 
in Sec.UV] Second, it is clear that (cos(9T)) n approaches 
exponentially with n for TO ^ pn, which is usually the case, 
as 9 is unknown. As a consequence, the SNR of the outcome 
is weakened, especially for n ^> 1. 

In this paper we propose a different method to perform 
multi-parameter estimation that achieves the QML scaling T, 
Interestingly, our method focuses on the evolution of observ- 
ables such as tensor products of Pauli operators (i.e., we work 



in the Heisenberg picture), rather than the evolution of the 
state of the probe itself. For this reason, in Sec. Ill Al we start by 
giving a brief description of Hamiltonian evolution in terms of 
Lie algebras. We then present an adaptive Bayesian estimation 
method to estimate single parameters at the QML with DQC1 
(Sec. IIICI l. Moreover, in Sec. [Ill] we show that, by apply- 
ing dynamical-decoupling techniques and different Suzuki- 
Trotter approximations, multi-parameter estimation can also 
be performed with DQC1. Here, we deduce that when the 
amount of short-time evolutions is considered a resource, the 
QML is asymptotically reached in the order of the approxima- 
tion. In Sec.[IV]we discuss the particular example of reference 
frame alignment and show that to estimate the Euler angles 
dynamical-decoupling techniques are not required. In SecfVl 
we discuss the reasons why DQC1 allows us to reach the fun- 
damental quantum limit in some cases, even though this model 
is less powerful than standard quantum computation. Finally, 
we present the conclusions in Sec. [VT] 

In the following, we ignore the effects of decoherence in 
our quantum algorithms and we assume that all experimental 
parameters can be controlled with arbitrary precision. 



H. SINGLE-PARAMETER ESTIMATION 

When the n-qubit probe interaction can be described by a 
Hamiltonian H = 9Hq, single-parameter estimation aims to 
return an estimate 9 of the unknown 9 at the highest preci- 
sion possible, for some given amount of resources. For fixed 
n and || Hq I) , our main resource is determined by the to- 
tal evolution time under H. Let W(T) = e~ i6H ° T be the 
(unitary) evolution operator induced by H , during a time in- 
terval T. Clearly, if W(T) acts non-trivially on some op- 
erator O, information about 9 can be gained by computing 
h(9T) = tr [W^(T)OW(T)0]/2 n e C for different values 
of T. The form of h(9T) can be obtained through the repre- 
sentation theory of Lie algebras (Sec. Ill At . Contrary to the 
example given in Sec. [I] h(9T) can be estimated using DQC1 
without controlling the operation W(T) (Sec lIIBt . Assuming 
that the accuracy in the estimation of h(9T) remains constant 
regardless of T, an adaptive Bayesian estimator that returns 9 
at accuracy 0(1/T) (i.e., the QML) can be built in some cases 
of interest. Furthermore, to avoid signal loss due to possible 
large values of n, we choose the operator O such that h(9T) 
does not depend on n. These points are studied and explained 
in more detail below. 



A. The Heisenberg picture 

If X = is an observable acting on n qubits, we define 
(Heisenberg picture) 



X(T) = W f (T)XW(T) , 



(5) 



with W(T) = e l9H « T as above. Assume now that we are 
given two observables Hi and such that, together with 
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Hq, they generate an su(2) Lie algebra. We obtain (see ap- 
pendixlAl 

Hi(T) = VK f (r)iJiVK(r) = cos(26»T)iJi - sm(29T)H 2 . 

(6) 

] = 

(7) 
(8) 

are 



If the operators are Schmidt pseudo-orthogonal (tr [HiHj S 
dSij ) we have 

cos(2(9T) = tr[Hi(T)Hi]/d , 
sin(2(9T) = -tr [H 1 (T)H 2 }/d . 

With no loss of generality we expand Hj 

Ej=i e^cr^, where a„ d = a 1 ^ ® • • • ® <r£ 
tensor products of Pauli operators (henceforth simply called 
Pauli products), - € {1, a^, o"j,, <r z }, and € K. are 
known coefficients ifSOTI . We obtain 

cos(20T) = e^'V'Hr [W^ (T)a^.iW(T)<r^ yl ]/d , 

(9) 

sin(20T) = - ^ e' 1 ' V< 2 tr [^(T)^,!^^)^^]/^ . 

(10) 

That is, cos(2#T) and sin(2#T) can be estimated at a fixed 
precision by L 2 runs of the circuit of Fig. Q] Each run returns 
an estimate of tr [U^j-^ ,j']/2 n , with unitary U^j-^'j' = 
(T)a Pi jW(T)a fJi ^ji . In Sec. lIICI we show how to estimate 
9 at the QML from the estimation of cos(26»T) or sin(20T), 
for different values of T. Although we are mainly interested in 
the fundamental scaling achieved by increasing T, it is worth 
noting that the scaling with L 2 can be largely reduced in some 
cases of interest (see Appendix lAl. 



B. DQC1 circuits 

We now show how to avoid the need of controlled C W{T) 
operations, when estimating parameters with DQC1. This 
is of great importance since C W{T) may not be available 
for our use. To show this, we focus on the estimation of 
tr [W^(T)a,j, )-7 -W(T)o- M / ,_,•/], as required by Eqs. © and ([TOJ. 
The circuit that accomplishes this task is shown in Fig. |2] 
Here, the C U operation of Fig. Q] has been replaced by 
C W^ (T) C (T MJ - C W (T)°a y ji , where each operation is con- 
trolled by the ancilla a. Nevertheless, note that one can ac- 
complish the same task even if the action of the operators 
W{T) and W^{T) is not controlled |21]. That is, 

(11) 

which is clearly the identity operator when a is in |0} a . In 
Fig. [3] we show a simplified circuit that allows us to compute 
the above trace. The last operation W'(T) is not included 
as it does not alter the measurement outcome. Also, W^(T) 
may not be an available resource. Thus, the circuit can be im- 
plemented using (known) elementary gates and the available 
time-evolution operator W(T) only. 
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in this case. Since every unitary is controlled by a (filled circles), 
including the evolution operator W(T), the execution of this 
algorithm may be unfeasible due to the nature of the problem. 
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FIG. 3: Simplified version of the DQC1 circuit of Fig. [2] The opera- 
tors C W(T) and W^(T) are avoided. The estimation of cos(29T) or 
sin(26T) [Eqs. l[9j and d lOb l is performed using available resources 
only. 



The next step is to detail a strategy to estimate 6 from 
cos(26*T) or sin(2#T), and to characterize the associated er- 
ror and resources needed. For this reason, we first make an 
assumption on the standard deviation A of the output returned 
by the DQC1 circuit of Fig. [3] Obviously, the smaller A is, 
the better the precision of the resulting estimation. We con- 
sider A <C 1, which can always be achieved by simple repeti- 
tion of the computation. In some cases, like liquid-state NMR 
quantum computation, where a vast amount of molecules con- 
tribute to the output signal, a A < 1 could be achieved 
in a single run. Our estimation procedure should take ad- 
vantage of this property by going beyond just performing a 
bit by bit estimation of 9, as it is done in several pure-state 
phase estimation techniques that involve strong (projective) 
measurements JH 0] . Moreover, we assume that A remains 
approximately constant in a certain region of values of 29T. 
This is a consequence of weak measurements l22ll . In view 
of the central limit theorem, we assume that the measure- 
ment outcome is distributed according to Af(cos(20T), A 2 ) 
or A/"(sin(20T), A 2 ) (see Sec. Hand Ref. |fl7tl ). 



C. Adaptive Bayesian estimation 

In Bayesian estimation, a parameter a to be estimated is 
considered to be a variable with an associated (known) proba- 
bility distribution f(a) (i.e., the prior distribution). The prior 
distribution formalizes the experimenter's state of belief about 
a. It is the job of the experimenter to gain access to a sample 
of data {xi, . . . , xk] whose distribution f(xjc, ■ ■■ ,xx\a) de- 
pends on a. Thus, the joint distribution of {x%, . . . , xk } and 
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a is 



/(%, ■ ■ ■ , xi, a) = f(x K , x 1 \a)f{a) . (12) 

After observing a set of measurement outcomes 
{xi, . . . ,xk}, this information is used to obtain a pos- 
terior distribution f(a\xx, ■ ■ ■ , x i) that corresponds to the 
experimenter's updated state of belief about the unknown a. 
This update is done using Bayes' rule 

f(a \x K , ■ ■ ■ , xi) = — 7} r- , (13) 

f{x K , ■ ■ • 

where the marginal sampling distribution f(xx, ■ ■ ■ , x\) can 
be calculated from the joint distribution as 



!{xk, ...,xx)= / f(x K , ■ ■ .,xi,a)da 



(14) 



In standard Bayesian estimation the probability of ob- 
serving an i.i.d. sample {x\,...,xk} is determined 
by the total sampling distribution f(xK,---,xi\a) = 
f(xj(\ce) • • • f(x\\a). In adaptive Bayesian estimation the 
outcome of the first measurement X\ can be used to con- 
trol the sampling distribution of the second measurement, 
f(x2\xx, a) (see below). In general, the sampling dis- 
tribution of the Zth measurement outcome can be condi- 
tioned by {xi, At the end, f(xx, ■ ■ ■ , xx\a) = 
f(x K \x K -i, . . . ,xt, a) ■ ■ ■ f(xi\a). 

The Bayes' risk quantifies the expected penalty to be paid 
when using a particular estimator &(xk, ■ ■ ■ , x\) of a, and a 
given cost function. It is common to search for an a that min- 
imizes the Bayes' risk with a quadratic cost function, given 
by 



(a - a(x K , . . . , xi)) 2 f(a\x K , 



, x\)da 



(15) 



This risk is just the variance of the estimator and is minimized 
by the expectation value of the posterior distribution, giving 
the optimal estimator 



&{x K ,---,x\) = I af(a\xK,---,x\)da. 



(16) 



While a more detailed explanation of the adaptive Bayesian 
procedure is given in AppendixlBl in the following we present 
a generic step of the estimation. Denote by A <C 1 the output 
precision of DQC1 when measuring cos(20T), which may ac- 
tually involve many (L 2 > 1) different runs of the circuit of 
Fig. [3] Assume that, from I previous estimations, we have ob- 
tained an estimator 9i of 9 such that, for known evolution time 
Ti and integer pi, 



29m » tt/2 + 2p i7 r 



(17) 



Furthermore, assume that the 95% confidence interval for the 
estimation is 

9i - 1.96A,/(2Ti) < 6 < §i + 1.96A,/(2Ti) , (18) 



with A; < A. Next, we show how to zoom in on 9i to obtain 
an estimator 6*;+i, so Eqs. (JT7J and ( fT8l are still satisfied when 
replacing I — > I + 1. 

To do this, we first find 7} + i such that 



20jT I+ i = k/ 2 + 2 Pi+^ » 



(19) 



with pi + i > pi integer. The (I + l)th measurement returns 
xi + i, an estimate of cos(a; +1 ), with ai + i = WTi +1 . This 
is done by running the algorithm of Fig. [3] with T = TJ +1 . 
We chose pi + i close enough to pi such that, for a' [+1 = 
— (tt/2 + 2p; + i7r), we approximate cosa/+i « — a' l+l 
(see AppendixlBl. Using Bayes rule and the joint distribution 
for the (I + l)th measurement, we obtain 

f(xi +1 ,ai +1 \xi,...,x 1 ) (20) 
= f{xi +1 \a l+li xi, ■ ■ ■ ,x 1 )f{ai +1 \x h ■■■ ,x{) 

(»i+i+"i+i) 2 



1 



1 



fry' I 2 



2ttA 



2nai +1 Ai 



with ai + i = T;_|_i/T;. Equation < f20b determines 
the posterior distribution /(aj+i xi, . . . , x\) = 
f(xt+i,at+i\xi,...,xi)/f(xi + i\xi,--- ,xi). This adap- 
tive procedure returns a new estimator f?;+i, and standard 
deviation A; + i/(2T/ + i), with Aj+i < A, satisfying 

9 l+1 - 1.96A, + i/(2T I+ i) < 9 < 9 l+1 + 1.96A, +1 /(2T I+1 ) . 

(21) 

The total number K of estimations is chosen such that the 
final standard deviation is reduced below the desired preci- 
sion. A sufficient condition is A/ (2Tk) < Ag. The fact that 
the standard deviation is reduced by Ti at each step [Eq. < f2TT >1 
guarantees that the QML is achieved (AppendixlBl. 

Remarkably, the confidence level of estimating the mean 
with error e in our algorithm increases exponentially as 1 — 
e -c(t/r) ) ( w j tn q > g an( j T ^ correS p 0n( ji n g standard de- 
viation. This is clearly an advantage with respect to the stan- 
dard pure-state phase estimation algorithm, where the confi- 
dence increases as 1 — 0(1/ e) lfl2ll . 



D. Black-box estimation: discrete time evolution 

Imagine now that, instead of being able to evolve under the 
action of H for any period of time T, we are given a black 
box whose action is to perform the unitary operation Wb = 
e~ lH . Like in the previous case, H = 9Hq- That is, we 
are only allowed to evolve under H for a discrete time by 
simple concatenation of Wg's. This condition restricts the set 
of accessible operations to elementary gates and operations of 
the form 



q times 
W(q) = Wb~^~W^ , 



(22) 



only. We seek to estimate 9 at the QML using a modification 
of the previous adaptive Bayesian method. 
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We now give the generic step for achieving the QML in the 
discrete time case (see Appendix ICl for the first step). We as- 
sume that 9i is the mean of the estimator obtained after the kh 
measurement, performed with qi E N* uses of Wb- Because 
we seek to make estimations around tt/2, we take 



1B x _ x q x + 24>i = tt/2 + 2 Pl ir, VZ 



(23) 



with 61-1 the mean of the estimation in the (I — l)th mea- 
surement, pi E N*, and <j>i the phase compensation. In 
other words, the kh estimation was performed by implement- 
ing the algorithm of Fig. with W(T) -> W a {qu4>i) = 
(WB) q 'e-^ lHo . Because we can always consider tt/2 > 
<j>l > —tt/2, the phase compensation is made with unit cost. 
We also assume that < A is the standard deviation of the 
estimator of 20qi Q . 

For the (I + l)th measurement, we write qi + i = bqi with 
small enough b 12411 . The (I + l)th measurement returns yi+i, 
an estimate of cos(/3;+i), with /3i+i = 29qi + i + 20;+i m 
tt/2 + 2np l+1 . Thus, to obtain f (j3i +1 \yi+i, ■ . .,yi) in the 
adaptive Bayesian step, we approximate 



f(m+i\yi, ■ ■ -,yi,0i+i) 



1 



2ttA 



(Vi+l+gj + l) 2 

-e 



with 



f3' l+1 = A+i - (tt/2 + 2 Pi+1 tt) . 
Moreover, since 

1 (<3 '+i 



(24) 



(25) 



■ • ■ ,yi) 



/2^6S, 



e 2 " 2 < s i> 2 , (26) 



the resulting distribution f(@i+i \yi+i , ■ ■ ■ , yi) is normal. This 
is a consequence of Bayes' rule. Its mean and standard de- 
viation, obtained by combining the exponents appearing in 
Eqs. (GUi and (|26]l, are 



*, {b' l+1 ? 

~i + (b> l+1 r yi+l ' 



= [b' l+1 /Jl + (b' l+i y A<A, 



(27) 



(28) 



with b' l+1 = 6S(/A < b. Equations d27l i and ( |28] > guarantee 
the success of the induction method. These quantities deter- 
mine the mean and standard deviation of the new (improved) 
estimator of 9 as 



V' 



1 

2b' 



tt/2 + 2pi +1 ir 
A 

< 



(&; + i) 2 



.-2/z+i 



2b l 2b l 



(29) 



The similarity of these results and those obtained for the con- 
tinuous time case [Eqs. ( IB20b and (1B21I) 1 is clear. 

Summarizing, if Ag denotes the desired precision in the 
parameter estimation, a sufficient condition is to choose the 
total number of measurements K such that the final precision 
satisfies A/(26 K ~ 1 ) < A e . Since the total number of uses 
of W B is given by 1 + • • • + b K ~ 1 = (b K - \)/{b - 1) = 
C[A/(2Ae)], the QML is also reached in this case. 



III. MULTI-PARAMETER ESTIMATION 

In this section we consider a more general case where the 
unknown interaction with the ?i-qubit probe can be described 
by a Hamiltonian 



H 



(30) 



Here, 9 V £ R and a v are Pauli products. Using the results 
of Sec. Ill CI we seek to estimate every (unknown) 9 V such 
that the returned precision approaches the QML for a given 
amount of resources or evolution time. For this reason, we 
assume that a previous estimate of every parameter, with mean 
tt > 9q > 0, is known. 

Using dynamical-decoupling techniques j2~5ll . the multi- 
parameter estimation case can be converted into P single- 
parameter estimations. For simplicity, we consider first the 
case where there is a Pauli product a such that (for some v) 



Then, 



[<r„,a] = , 
{£T„/,£r} = OVv'^v 



R v = 9 v a v = (H + aHa)/2 . 



(31) 
(32) 

(33) 



(In general, methods like the one just described can be used to 
decouple any H v from any H.) We define S V (T) = e~~ lHuT 
to be the corresponding evolution operator. If such an operator 
were to be an available resource, 9 V could be estimated using 
the scheme of Sec. Ill CI To do this, we would have to replace 
W(T) -> S U (T) and cr MJ -> a x in the circuit of Fig. [3] with 
{ov, <7i} being Pauli products that generate an su(2) algebra. 

We now show how to approximate S U (T) from accessi- 
ble operations that include W(T) and elementary gates only. 
To show this, we use a Suzuki-Trotter approximation 12611 . 
Specifically, for q = 1/e E N* , we decompose 

q times 



S v {T)=S v (eT)---S 1/ (eT) . 

If S„(eT) denotes a pth order Suzuki-Trotter approximation 
to S l/ (eT), we have 



with || 
Then, 



|| S u (eT) - S u (eT) ||= 0[\\ H f (eTf] , (34) 
|| some operator norm (e.g., the largest eigenvalue). 



e/2 =\\ S V (T) - S V (T) ||= 0[\\ H \\ p e^T?] , (35) 

with S V {T) = [S v (eT)] q . Equation d35j was obtained using 
Eq.®, together with || || = || (A-B)A^- 1 +B(A- 

B)A4~ 2 + ... + Bi- 1 {A- B) || < q || A - B ||, for A and 
B unitaries. 

We now show how to build S V (T) out of available re- 
sources. In the simplest case (i.e., p = 2) the evolution op- 
erator at short times factorizes as 



§ v {eT) 



e -iHeT/2 ae -iHeT/2 a 



(36) 
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Then, S y {T) can be implemented using H evolutions and a 
gates only. Similarly, a second-order Suzuki-Trotter approxi- 
mation (i.e., p — 3) is given by 

S u (eT) - e- iH ^' A cje- iH ^' 2 ae- iHcT / A . (37) 

Higher order approximations can be constructed in a similar 
fashion ll27ll . so they can always be implemented with accessi- 
ble gates. The larger p is, the shorter the actions of H in each 
step. Thus, we require precise time control in our algorithms 
(measurements). 

Replacing W(T) — > S V (T) and a^.j — ► <j\ in the circuit 
of Fig. fallows us to estimate 9 V ' . To show this, consider the 
measurement output z% obtained in the Zth measurement. Zi 
will give us an estimator of the angle 6 V plus a correction 

cos(2eT,)+7(e,J»,r I ,"?) > (38) 

with ~6 = 9 1 , ...,9 P . Thenormof7(e,p,T i ,"0 ) E Mean be 
bounded above as 

|7(e,p,7],^)| = 

2" 

^ tr || St{Ti)aiS v {Ti)ai - Sl(Ti)<nS v (Ti)ai \\ ^ 
— 2 n — 

(39) 

We have used Eq. d35l ) to obtain Eq. (l39l . Because 

j(e,p, Ti, 9) depends on the fl^'s, we consider it a vari- 
able with an associated (worst-case scenario) prior distribu- 
tion given by /( 7 ) = W(0, A 7 ), with A 7 = O(Ae) Q. 
The net effect in the adaptive Bayesian procedure is that now, 
the joint distribution determines a marginal distribution after 

7(e,p, Ti, 6) is integrated out. More precisely, for the Zth es- 
timation, we have 

f(xi\xi-i, . . .,xi,a>i) 

I f(xt\xi-i, . . . ,xi,ai,^)f(-/)d-y , (40) 

with a" = 29 v Ti. Making a linear approximation in the co- 
sine function, we obtain 

f{xi\xi- u . . . ,11, aH « M{a\ v ,(A') 2 ) , (41) 

with oil = a i _ (^/^ + Zpin) (see Sec. Ill Cb . and updated 
variance 

(A') 2 = A 2 + A^ . (42) 

Thus, we can use the adaptive method of Sec. lIICI to estimate 
every parameter 9 V at the QML, by replacing A — > A'. 
Notice that one could also implement an adaptive Bayesian 

approach to learn about j(e,p, Tj, 9 ). In such a case, the 
distribution f(-f) would need to be updated after each mea- 
surement based on the previous measurement outcomes. Nev- 
ertheless, we did not consider this approach in the above dis- 
cussion because we assumed that the Suzuki-Trotter approxi- 
mation used is good enough for our purposes (i.e., A 7 <C 1). 



When the resource of interest is the number of calls to 
S v (eT), the amount of resources to reach a precision Agv 
changes. For this reason, consider the total evolution time 
T t = 0[(A + A 7 )/A e ,)]. with A 7 = 0(A \\ H f e^Tf). 
Assume that we want to keep T t constant, regardless of p. 
Then, q=l/e = C[(T t ) (p/p_1) ]- That is, C[l/(A e ,) (p/p_1) ] 
actions of S v {eT) are required to attain Ag», and the QML is 
asymptotically reached in p. 

IV. REFERENCE FRAME ALIGNMENT 

Imagine that two distant parties, Alice and Bob, suffer some 
frame misalignment, which is manifest in the way they char- 
acterize their operations on equivalent quantum systems. This 
might be a result of not sharing synchronized clocks or hav- 
ing different spatial reference frames. Aligning both frames 
requires the exchange of physical systems carrying "unspeak- 
able" information. This information is encoded as frame- 
dependent parameters that need to be estimated [29]. The 
resource that limits the quality of this estimation is the num- 
ber of systems interchanged between Alice and Bob. We will 
show that Alice and Bob can align their frames within the 
DQC1 model at the QML. In particular, we propose a modifi- 
cation of a pure-state protocol for frame synchronization 
based on repeated coherent exchanges of the n-qubit probe 
only. Remarkably, the state of the probe remains completely 
mixed and separable from the ancilla at every step Qj| . More- 
over, in our protocol, Bob never accesses the ancilla that Alice 
measures. 

We consider first the case where the effect of the frame 
misalignment is uni-parametric. That is, Alice's and Bob's 
description of operators acting on equivalent Hilbert spaces 
is known to differ by a unitary transformation Vg = e °, 
with 9 unknown. More explicitly, for some operator O we 
have 

O b = VfjOVe = e ieH °Oe- ieHo , (43) 

where we used the superscript B to denote the action of O 
in Bob's frame. Since Hq is known, we assume that we can 
find pseudo-orthogonal observables Hi and H2 such that they 
form an su(2) algebra [i.e., Eq. JA5b is satisfied in Alice's 
frame] . 

An elementary step of the protocol starts with Alice sending 
the n-qubit probe to Bob. Subsequently, Bob applies the oper- 
ation e~™ Hl / 2 and returns the probe to Alice. Finally, Alice 
applies the adjoint operation e l7lH i/ 2 . The resulting operation 
on the state of the probe in Alice's frame is 

e mH 1 /2 e -iwH?/2 _ /2 e i9H ^iitHt /2 e -i9H Q 

= e - 2ieH ° = V-28 ■ (44) 

Equation ( PPfl i can be obtained by working in any faithful rep- 
resentation of su(2). The global action of each step can be 
seen as a "black box" implementation of V20, whose param- 
eter we want to estimate. Using the results of Sec. Ill Dl the 
circuit of Fig. [3] can be used to make a first estimate of 9 if 
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W{T) is replaced by V 2 e (see Fig. To zoom in on previous 
estimations requires instead the implementation of the unitary 
V 2m g = (V28)" 1 , with m G N*. This can be done by sim- 
ple concatenation of elementary steps, requiring m coherent 
exchanges of the probe. 



0> a a <0| ^ 

P 



ln/2 



Bob 



FIG. 4: Elementary step of the frame alignment protocol. The n- 
qubit probe, whose state remains completely mixed, is exchanged 
between Alice and Bob. The Pauli products cr M j and cr M / ji appear 
in the decomposition of Hi [see Eqs ([9} and dlOH. 

Another case of interest is the alignment of spatial reference 
frames |5]. Let us assume that Bob's operators are related to 
Alice's through a rotation 



R _ e -iil>H 2e -ieH le -i<t>H 2 



(45) 



where {</>, 8, tp} are Euler angles and {ifi,/^} generate an 
su(2) algebra. Consider a synchronization protocol with this 
elementary step: Alice sends the probe to Bob; Bob applies 
the operation e~™ H ^ / 2 and returns the probe to Alice; Alice 
applies the operation e mH2 / 2 . The effective rotation of this 
step, in Alice's frame, is given by 



V 



rff 2 /2 i jt e -wff 2 /2 jR . (46) 



The action of V on H 2 implies that 

tr [V' f H 2 V'H 2 ] = dcos(26>) , (47) 

with d some normalization constant [Eq. (9)]. Again, the pre- 
vious derivation can be carried out in any faithful represen- 
tation of su(2). Thus, a first estimation of the Euler angle 
can be performed by using the circuit of Fig. [3j replacing 
W(T) — > V', and with j being the Pauli products appear- 
ing in the expansion of H 2 . Furthermore, since 



tr [{V*) m H 2 (V') m H 2 ] = dcos(2r 



(48) 



we can zoom in on the previous estimation by applying V', 
m > 1 times, and using the results of Sec. Ill Dl In this case, 
m coherent transports of the probe between Alice and Bob are 
required. 

Notice that in this estimation procedure, the action of rele- 
vant operations need not be controlled by the ancilla. Finally, 
the other Euler angles can be estimated in a similar way if 
Alice and Bob agree to apply other analogous operations. 



V. PARAMETER ESTIMATION VS. GROVER'S SEARCH 
ALGORITHM IN DQC1 

So far we have shown that to reach certain precision in the 
estimation of a parameter, DQC1 requires less resources than 



other methods. One may wonder if such a quantum speed- 
up can also be attained in problems such as searching for a 
particular property in a given set (i.e., search problem), which 
is the case for pure-state quantum algorithms. However, in 
Ref. lITill the authors proved that DQC1 is strictly less power- 
ful than standard quantum computation in the oracle setting. 
This implies that DQC1 does not provide a quadratic quantum 
speed-up in the search problem as the one given by Grover's 
algorithm j30h . To show this, consider the situation in which 
we are given a black box that implements either the unitary 
U B = e ie \ s ^ s K with 6 ^ 0, or U B = X over the state of 
the n qubits in the probe. Here, IS") encodes the solution to 
our search problem. We want to determine the existence of 
a solution; that is, we want to specify if the action of Ub is 
trivial (i.e., 11) or not (i.e., e l6, ' s ^ s '). In fact, this is a phase 
estimation problem in which we have to decide whether the 
phase is or 6. 

In general, the output of a DQC 1 algorithm is given by 



tr [pftrl 



(49) 



with p f the ancilla-probe state right before the measurement. 
With no loss of generality, 



Pi 



WqU b ■ ■ ■ WxUbWoPqWIuIwI ■ ■ ■ U B WX , (50) 



with po = (|0) a a (0| ® ll„)/2 n being the initial state, and Q the 
number of calls to Ub- Since |0) a a (0| = (l a — o"|)/2, we have 

(o* g ) = tr [WqU b ■ ■ ■ U B W alWlU B ■ ■ ■ U B W^Qa 3 z }/2 n+1 

(51) 

Following the proof in Ref. II 1411 . we obtain 



{<)\u B =t- {al)\ UB=e is\ S){S \\ <4Q/2 



TI + l 



(52) 



Since DQC1 returns (<r a z ) at accuracy A, Q must be exponen- 
tially large in n or the algorithm needs to be executed expo- 
nentially many times to determine whether there is a solution 
or not. That is, if J is the amount of times that the algorithm 
is executed, we would expect that the precision in the estima- 
tion scales as A' = C(A/\/j). To solve the problem, it is 
necessary (but not sufficient) to choose J and Q such that 

4Q/2 ,l+1 > A' , (53) 

requiring = JQ uses of the black box. Thus, 

N > 0(2") . (54) 

If Q = C(\/2") as in Grover's algorithm, we need J = 
0(\/2™) to satisfy Eq. d54l ), and no quantum speed-up with 
respect to the classical counterpart is obtained in this case. 

The reason for the existence of a quantum-speed up in pa- 
rameter estimation is that the unitary operators considered 
act non-trivially in a large-dimensional subspace of the corre- 
sponding 2™-dimensional Hilbert space. In this case, the out- 
put signal obtained by executing the DQC1 circuits enables us 
to distinguish between different unitaries and to estimate the 
parameter. In the search problem, however, the operator Ub 
is very close to the identity operator 11 in that it only affects 
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the state \S). Therefore, its action over highly-mixed states 
is almost trivial. This is not the case in pure-state algorithms 
where one usually works in a two-dimensional Hilbert space 
spanned by the states \S) and jS"- 1 -), with (S\ S- 1 ) = 



VI. CONCLUSIONS 

Parameter estimation at the quantum metrolo gy l imit could 
have a wide range of applications in metrology llOl [Till . Fur- 
ther, single-parameter estimation is related to phase estima- 
tion, a cornerstone of quantum computation. Mixed-state 
quantum computation, as formalized in the DQC1 model, is 
interestin g bo th from a theoretical and from a practical point 
view Il4l ll5[ [Till . We have shown that, under fairly general 
conditions, it is possible to perform parameter estimation at 
the quantum metrology limit within the DQC1 model. These 
conditions are presented using Lie algebraic methods. The 
algorithm proceeds using adaptive Bayesian estimation. In 
each step we zoom in on the previously estimated parameter, 
while ensuring that the increased variance remains below cer- 
tain bounds. A measurement reduces such a variance to the 
previous value and this procedure is repeated until the desired 
precision in the estimation is reached. In sort, the procedure 
ensures that the phase is kept in the same region, with almost 
constant variance, but increasing winding number. 

The adaptive estimation is clearer when the time of the 
evolution under the unknown Hamiltonian (parameter) can be 
controlled at will. When lacking this freedom, as in the case 
of black box estimation, the algorithm for continuous time can 
be amended with some straightforward modifications. More- 
over, to perform multi-parameter estimation, implementation 
of dynamical-decoupling techniques reduce the problem to 
several single-parameter estimation procedures. Yet, these 
techniques are not necessary in some simple cases like frame 
alignment between two parties, Alice and Bob. For spatial 
frame alignment, the Euler angles can be estimated at the 
quantum metrology limit with Bob having access to the com- 
pletely mixed, separable, state of the probe only. Surprisingly, 
although precise estimation is intimately related to the quan- 
tum speed-up given by Graver's search algorithm, the later 
cannot be performed efficiently within the DQC1 model. 
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APPENDIX A: LIE ALGEBRA REPRESENTATIONS 

For X(T) = W^(T)XW(T) we can also write X(T) = 
X + iT[H,X]-T 2 /2[H, [H,X]] + - • -,with [Y,X] = YX- 
XY. That is, X(T) is a linear combination of observables 
belonging to the Lie algebra rj generated by H and X lfl9tl . 
In general, f) = {0\, . . . , Om} is an -A/-dimensional (real) 
semisimple Lie algebra, with Oj = Cr, and tr [OjOj] = dSij 
(d £ R). A faithful representation of f) is the mapping Oj — > 
Oj, with Oj = Oj being (m x to)— dimensional matrices that 
satisfy 

Pi, Oj] = £ f^k ^ [Qi, Oj] = £ f:\(K , (Al) 
fc fc 
tr [O t O j] = dS l3 , (A2) 

for some d £ R. The coefficients are the so-called struc- 
ture constants of h. Then, X(T) = J2 k c k O k , c k £ R, and 

X(T) = W\T)XW{T) = c k O k , (A3) 

k 

with W(T) = e- lRT . The c fe 's depend only on the f k, s. 
Equation dA3t implies that if two different sets of (linearly- 
independent) matrices have the same commutation relations, 
the coefficients c k can be determined by working in either ma- 
trix representation: 

c k = tr [X(T)O k ]/d = tr [X{T)O k ]/d . (A4) 

Because {Hq, Hi, H2] span a su(2) Lie algebra, they sat- 
isfy 

[Hj,H k ]=2ie jkl Hi, (A5) 

with j,k,l £ {0,1,2} and Cj k i the totally antisymmetric 
symbol. The su(2) Lie algebra can be built upon (2 x 
2)— dimensional Hermitian, traceless, matrices (i.e., Pauli 
spin- 1/2 operators). This allow us to carry out the calculation 
of Eq. iA3i in a low-dimensional representation. 

To give an example where the L 2 trace estimation can be 
reduced [Eqs. d9|10H , consider again the situation where H = 
2j=i °z [Eq- ©]. Then, for example, 

al(T) = e iHT ale- iHT = cos(26T)<j 1 x - sin(20T)c^ , 

(A6) 

as [a^ai] = V j > 1. That is, cos(26>T) can be es- 
timated by computing the renormalized trace of the unitary 
U = cr^(T)cr^. only. This situation can be generalized to the 

case when Hq = Ylu=i e>1 ' 0(7 nfl> if there is an Hi = a x , with 
a x a Pauli product such that (for some fi) 

[0^,0, <V,o] = V/i,/i', 
[<V,o,<^] = V/i'^M- (A7) 

It follows that 

Wij.,0, Ci] = 2icr 2 , 

[<V,o.*J = V/i'^, (A8) 
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with <r 2 = —i(<J ll aa , i). Equation (IA8b results from the Jacobi 
identity and the definition of er 2 . Thus, {(Tuft, c 2 } satisfies 
the 5u(2)-commutation relations, yielding 

cos(20e M '°T) = tr[W t (T)a 1 W(T)a 1 ]/2 n , (A9) 
sin(20e M '°T) = tr [M /t (T)CT 1 VF(T) ( T 2 ]/2" . (A10) 



We estimate Eqs. (|A9t and (lAlOb by executing the cir- 
cuit of Fig. Q] with t/ = W^(T)a 1 W(T)a 1 and [/ = 
W^(T)a±W(T)a z , respectively. 



APPENDIX B: THE ADAPTIVE BAYESIAN ESTIMATION 
PROCEDURE 

To reach the QML in our estimation procedure, we first as- 
sume an initial estimation of 9 given by a prior distribution 
AT(0q, A'). For simplicity, in the following we focus on the 
case where 9 is determined from the estimation of cos(20T) 
for different values of T, denoted as TJ. This is done using 
the DQC1 algorithm of Fig. [3] where the corresponding Pauli 
products er„ j are determined by Eq. (|9). Thus, X[ denotes an 
estimate of cos(20Ti) and we consider n > 9q > 0. Other- 
wise, the estimation of sin(20Ti) is also required to determine, 
for example, the corresponding quadrant of 9. 

To obtain xi, we choose Ti such that 20 o Ti = n/2 l3lll . 
This can be done with an upper bounded initial use of re- 
sources when O € [x, t), with \ > 0. That is, 1/4 < T\ < 
7r/(4x). Nevertheless, if 9q <C 1 (i.e., Xi ^> 1), a similar anal- 
ysis as the one carried out below can be performed by mea- 
suring sin(20T{) instead of cos(20Xi), with T[ = 0(1) jH. 
Therefore, we take Af(tv/2, (c A) 2 ) as the prior distribution 
of «i = 29Ti. We define here A to be the output precision 
of DQC1 when measuring cos(20T), which may actually in- 
volve many (L 2 > 1) different runs of the circuit of Fig. [3] 
We assume 1 ^> cA A. The measurement outcome x\ of 
the first measurement has then a sampling distribution given 

byA/"(cos(ai),A 2 ). 

The joint distribution /(.Ti, «i) = /(xi|ai)/(ai) is 



f(xi,ai) 



1 



i'i-™i°i)f 



'2ttA 
1 

^a 



2ttcA 



o»i+°>ir 



1 



K) 2 

■;,:2 \2 



2?rcA 



(Bl) 



Here, a[ = ot\ —n/2 and, for simplicity, we approximated 
cos(ai) at first order by — a[ so that the joint distribution 
is normal. The error in this approximation can be bounded 
above, with high confidence, as 



| cos(ai) - (-ai)| < S = (c'A) 3 /6 « 1 



(B2) 



for some c' > c. For example, if choosing d = 1.96c, 
Eq. dB2b determines a 95% credible interval for cos(«i). Such 
a confidence can be made exponentially close to 1 as c! in- 
creases. Of course this error can be avoided if no approx- 
imation is made and other analytical or numerical methods 
are used. Nevertheless, a linear approximation to the cosine 



is enough for our purposes, as it will yield the proper re- 
sults Il33ll . Moreover, the error of the above approximation 
will be further corrected by subsequent measurements. This 
is a consequence of the adaptive method. 

The following step is to update the information about ol\ 
(or 9) based on the measurement outcome x\. The posterior 
distribution is /(ai|a;i) = f(x\,a.\)/f(x{). Using Eq. (IB lb 
this distribution can be shown to be 

/(ai|xi) «AA(«i,A 2 ) . (B3) 

&i and Ai are determined from the exponent E = —(xi + 
ai) 2 /(2A 2 ) - ai 2 /(2c 2 A 2 ) in Eq. (ED- We write 



E 



2A 2 



+ g(xi), (B4) 



implying 



Ai = 



VT 



:A < A , 



&! = ir/2 



l + c 



(B5) 
(B6) 



Summarizing, if 9\ is our estimator of 9 after the first mea- 
surement, we have 



cti 1 / c" 
n = — = — \ k/2 - xi 1 , 



2Ti 2Ti 



l + c 2 



(B7) 



which is only a linear correction in X\. Moreover, our knowl- 
edge about 9 has increased such that 

0i - 1.96Ai/(2Ti) < 9 < 0i + 1.96Ai/(2Ti) (B8) 

is a 95% credible interval. If the desired output precision Ag 
in the single-parameter estimation satisfies Ai/(2Ti) < Ag, 
the estimation procedure stops here. Otherwise, further mea- 
surements are required. 

In the previous analysis we have neglected other values 
of ot\ mod 2tt that would yield the same measurement out- 
come. This assumption introduces an extremely small error 
bounded above by 



erfc 



2tt \ 



;A\/27 



(B9) 



where erfc denotes the complementary error function of 
the normal distribution. Since cA <C 1, we have 
erfc(27r/ (cA\/2)) ~ 0. For example, for the unrealistic case 
of cA = 1, we have erfc(v / 27r) ~ 3.4 10 -10 . Thus, our 
assumption does not affect the final results of the estimation 
procedure. 

To increase the precision we zoom in on the previously 
estimated 0, so a better estimate is attained. Since x\ = 
0{cA) < 1 JH, we have \26iTi - tt/2\ < 1. Accordingly, 
there exists a time period T2 > T\ such that 



20i T 2 = tt/2 + 2p 2 7r 



(BIO) 
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with p 2 G N*. T 2 denotes the evolution time in the second 
measurement (estimation). The main reason behind Eq. dBlOt 
is that here, as in the first measurement, we seek to make a 
phase estimation around tt/2 mod 2ir. In this region, the co- 
sine function is more sensitive to variations in the phase. The 
second measurement returns X2, an estimate of cos(a 2 ), with 
a 2 = 29T 2 . 

Since T 2 = a 2 T\, with a 2 > 1, we obtain a 2 ps (1 + 
4p 2 ) [see Eqs. ( IB7b and ( IBIOM . The previously estimated 8 
has a variance of order A/(2T 1 ) [Eqs. jB5i and dB8l)l, so the 
variance of a 2 is of order a 2 A. To guarantee that this variance 
is similar to that of ot\ (first measurement), we choose a 2 as 
large as possible so that IT35ll 



< d - 4 < a 2 < d 



(Bll) 



This implies that Eq. ( IB2b is still satisfied when replacing a.\ 
by a 2 , and a[ by a' 2 = a 2 — (7r/2 + 2p 2 7r). As an example, 
consider the case when d ps 10. Therefore, we choose a 2 ps 
9 < c', corresponding to p 2 = 2 in Eq. ( IB 1 Ob . 

After the measurement, the outcome x 2 is used to update 
our information about 6. Since d A « 1, 02 € [2p 2 7r, (2p 2 + 
1)7t] with large confidence. That is, we do not consider other 
values of a 2 mod 2ir and we make the estimation in this re- 
gion only. The joint distribution is now 



f(x 2 ,a 2 \xi) = f{x 2 \x ll a 2 )f(a 2 \x 1 ) 

c* 2+a ' 2 ) 2 1 — 



1 



2ttA 



1 

27ra 2 Ai 



(B12) 



f(a 2 \xi) has been determined using Eq. dB 10b . Thus, 
the posterior density distribution f(a 2 \x 2 ,xi) = 
f(x 2 ,a 2 \xi)/f(x 2 \xi) is a normal Af(a 2 , (A 2 ) 2 ), where 
d 2 and A 2 are determined by the exponent appearing in 
Eq. (IB 12b . These are 



Q'2 



A 2 



(tt/2 + 2 P2 TT 



1 + (4)2 

A < A , 



,r 2 . 



(B13) 
(B14) 



with 



a' 2 = a 2 Ai/A < a 2 . (B15) 
Summarizing, the estimator after the second measurement is 



? 2 = ^r( 



K) 



1 + K) 2 



•'2 



(Bi6) 



with a 95% credible interval 

e 2 - 1.96A 2 /(2T 2 ) < 6 < 6 2 + 1.96A 2 /(2T 2 ) . (B17) 

The standard deviation in the estimation of 9 has been reduced 
by a factor of order a 2 with respect to the one returned in the 
first measurement [Eq. flB8H . 

If the desired output precision satisfies A 2 /(2T 2 ) < Ag, 
the estimation stops here. Otherwise, we continue with the 
adaptive procedure. At each step I we find T; such that 



26i-iTi = tt/2 + 2 Pi tt 



(B18) 



where istheBayes' estimator determined by the previous 
measurement outcomes. Since pi > pi-\ are positive inte- 
gers, we write Ti = offi-i, with ai w (l + 4jjj)/(l + 4pj_i). 
The Ith measurement returns xi, an estimate of cos(a;), with 
ai = 28Ti. This is done by running the algorithm of Fig. [3] 
with T = Tj. To keep the variance of ai at order c'A we 
choose a; such that c' — 4 < a/ < d IT35ll . With this choice, 
Eq. ( IB2b is still satisfied when replacing ai by a/, and ctj by 
a[ = ai — (tt/2 + 2piit). Using Bayes' rule, and considering 



f(xi,ai\xi-!, 
1 



,xi) « 

(x !+ „;)2 



1 



/2^A 



2?ra;A 



we obtain for the Zth estimation 



20,T, = 7r/2 + 2p i7 r- 



(-1 



W) 2 



A;_ 

2T ; 



A 



K) 2 

A 



•'7 



< 



with 



y/\ + (a'J 2 m 2T ; 



a;A;/A < a ; 



(B19) 

(B20) 
(B21) 

(B22) 



That is, the variance of the /th estimator has been reduced by 
a factor l/Tj. 

We now show that the QML has been achieved. Consider 
the total evolution time T t = T\ + ■ ■ -+Tk, for K estimations, 
with T z = fllLaOiOTi. Then, 



Ti 



+ 



1 



T K ■ 



a 2 . . . a K 

Moreover, since (1/a;) < l/(c' — 4) < 1, we have 

4. 



T t < l(c' - 4)-( if - 1 ' 



l]Tjf < 



A" 



(B23) 



(B24) 



Equation (IB24b gives the total resource scaling T t = 
0[A/(2Ag)] (i.e., the total evolution time under the action 
of H), implying that the QML is attained. 



APPENDIX C: FIRST STEP IN BLACK-BOX ESTIMATION 

For simplicity, assume 7r/4 > #0 > 0, with 60 the mean 
of the prior (normal) distribution of 9. Then, there exists a <fii 
such that 



2§ + 20i = tt/2 



(CI) 



The first estimation of 6 can be done by running the algorithm 
of Fig. © replacing W{T) -> W o (l,0i) = %e-* fl ». 
Since i^o is known, W (l, <t>i) can be implemented with avail- 
able gates. The output of the first measurement, denoted by 
2/1, is a measurement of cos(/3i), with /3i = 2(0 + <f>{). 
The a priori distribution of fti is then given by f((3\) = 
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Af(ir/2, (cA) 2 ) ll36ll . Similar to the continuous time evolu- 
tion case, this distribution yields the joint distribution 



f(yi,Pi) = 



i 



1 

^A 



2A^ 



e 2A^ 



1 



V27TCA 

1 (/3l) 2 



2ttcA 



(C2) 



with f3'i = [}\ — tt/2. The first measurement returns an es- 
timator of Pi, with mean /3 X = tt/2 — [c 2 /(l + c 2 )]^!, and 
standard deviation Ei = (c/Vl + c 2 )A < A. Therefore, the 



first estimation of 9 has mean and variance determined by 

^1 = 5(^/2-^^)-^, (C3) 

Ei = Ex/2 = -7=^=^ < A/2 . (C4) 

Clearly, the accuracy of the estimation has increased after 
the first measurement. Thus, we continue with the adaptive 
Bayesian method by zooming in on the previously estimated 
parameters. 
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